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Abstract 



In some scenarios, the peculiar gravitational potential of linear and mildly 
nonlinear structures depends on time and, as a result of this dependence, a late 
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integrated Sachs- Wolfe effect appears. Here, an appropriate formalism is used 
which allows us to improve on the analysis of the spatial scales and locations of 

X 

the main cosmological inhomogeneities producing this effect. The study is per- 
formed in the framework of the currently preferred flat model with cosmological 
constant, and it is also developed in an open model for comparisons. Results 
from this analysis are used to discuss the contribution of Great Attractor-like 
objects, voids, and other structures to the CMB anisotropy. 

Subject headings: cosmic microwave background — cosmology:theory — large-scale 
structure of the universe — methods:numerical 
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1 INTRODUCTION 

In some scenarios where linear and mildly nonlinear structures create a time varying 
gravitational potential, the photons of the Cosmic Microwave Background (CMB) un- 
dergo a late Integrated Sachs- Wolfe (ISW) effect. In the absence of any cosmo logical 
constant, the partial time derivative of the gravitational potential tends to zero as the 
universe approaches a flat one and, consequently, the ISW effect tends also to zero. 
This paper is devoted to the study of some aspects poorly known of the ISW effect: 
our goal is a detailed analysis of the locations and scales of the subhorizon structures 
contributing to this effect. We are particularly interested in the scales corresponding 
to observable objects as voids, the Great wall, et cettera and, by this reason, we will 
only consider scales smaller than the horizon. We choose an adequate formalism to 
deal with this analysis. 

Two scenarios are considered: a flat universe with cold dark matter (CDM) and 
cosmological constant and an open universe with CDM, they are hereafter referred to 
as scenarios (or models) I and II, respectively. In both cases the spectrum has been 
first normalized by the condition as — 1, to consider other normalizations (other <jg 
values) when necessary. In case I, the spectrum corresponds to cold dark matter 
(CDM) with Q d = 0.25, tt b = 0.05, Q x = 0.7, h = 0.65 and n = 1, where Q b , tt d , and 
Q\ are the density parameters corresponding to baryonic matter, dark matter, and the 
cosmological constant, respectively, h is the reduced Hubble constant (h = H /W0, 
H being the Hubble constant in units of Km/ s.Mpc), and n is the spectral index 
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of the primordial scalar energy density fluctuations. The scenario II involves CDM 
and the relevant parameters are: Q d = 0.25, fi 6 = 0.05, h = 0.65 and n — 1. Model 
I is currently preferred according to recent observations of far la supernovae and 
the CMB spectrum (location of the Doppler peak), while model II can account for 
the abundances of rich clusters and Einstein's rings and, here, it is mainly used for 
comparisons. As it is well known, the normalization a$ = 1 does not lead to a good 
normalization of the CMB angular power spectrum in most cases; in other words, 
when the Ce coefficients are calculated (for <r 8 = 1), the resulting values do not fit 
well with the values observed by COBE, TENERIFE and other experiments. In each 
scenario, appropriate fits to the observed CMB spectrum correspond to er 8 values 
which are, in general, different from unity. The so-called bias parameter is b — l/cr$. 

Since our attention is focused on subhorizon scales, we will estimate the late 
ISW anisotropy in the ^-interval (10,40). For £ < 10, super-horizon scales would 
be important (see Kamionkowski and Spergel, 1994) and, then, the spatial curvature 
could be only neglected in model I; furthermore, the cosmic variance would lead to 
important uncertainties (ACe/Ce is proportional to [2/(2£+l)] 1//2 , see Knox, 1995) and 
the Sachs- Wolfe effect would be very important. For £ > 40, the Doppler effect starts 
its domination hidding other effects as the ISW one. An appropriate linear approach 
is used in next sections to estimate the Ce coefficients for 10 < £ < 40. The method 
used to do this estimation should facilitate the separation of the ISW effect from other 
contributions to the angular power spectrum and, moreover, this method should 
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give information about the sizes and locations of the main subhorizon structures 
contributing to the late ISW effect. The numerical integration of the Boltzmann 
equation or the computational strategy of Hu and Sugiyama (1994) could be used to 
perform the analysis of this paper; nevertheless, another appropriate approach -based 
on a certain approximation to the sources- is described and used in next sections. 

In previous papers, it was claimed that some Great Attractor-Like (GAL) objects 
located between redshifts 2 and 30 in open enough universes (without cosmological 
constant) could account for an important part of the Integrated Sachs- Wolfe (ISW) 
effect. Arguments in those papers were based on the Tolman-Bondi (TB) solution 
of Einstein's equations, which was used to estimate both the anisotropy produced 
by a single GAL structure (Arnau, Fullana & Saez, 1994; Saez, Arnau & Fullana, 
1995) and the abundance of these structures (Saez & Fullana 1999). Unfortunately, 
our TB simulations have some features, as the spherical symmetry and a particular 
form of compensation, which could affect abundance and anisotropy estimations. By 
this reason, the mentioned claim should be discussed using a general formalism (not 
TB solution). It is done in Section 4 as a subsidiary application (in model II) of the 
formalism described along the paper. 

In section 2, the method used to compute the angular power spectrum inside the 
t interval [10, 40] is described. Results are presented in Section 3 and, Section 4 
is a general discussion and a summary of conclusions. Finally, some words about 
notation: whatever quantity "A" may be, A L and A stand for the A values on the 
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last scattering surface and at present time, respectively. Simbols x\ 0, v, n, p B , 
S, a, t, G, stand for the comoving coordinates, the peculiar gravitational potential, 
the peculiar velocity, the unit vector in the observation direction, the background 
mass density, the density contrast, the scale factor, the cosmological time, and the 
gravitational constant, respectively. Units are chosen in such a way that the speed of 
light is C — 1. Quantities Vt Q and Vt m are defined as follows: Q = Q b + Q d + Q\ and 
Q m = Qb + ^d- The comoving wavenumber is k c , while k is the physical one. 

2 Ct ESTIMATIONS IN THE t INTERVAL [10,40] 

We are interested in the ISW effect produced by structures much smaller than the 
horizon scale and, consequently, the region of the hypersurfaces t = constant where 
the inhomogeneities interact with the CMB can be considered as flat; namely, the 
spatial curvature can be neglected in the open model II. This means that, even in the 
open case, the spatial part of the functions defining the linear structures under consid- 
eration can be expanded in plane waves. It is not necessary the use of the complicated 
solutions of the Helmholtz equation, which should be used in open backgrounds to do 
an exact and rigorous treatment of structure evolution. Of course, in model II, the 
time evolution is studied taking into account the existence of a space-time curvature 
distinguishing the open universe from the flat one. 

The potential approximation is used in our estimates. The basic equations (Martinez- 
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Gonzalez et al. 1994, Sanz et al. 1996) are: 

AT 1 r f ° 86 

— = - 4 , L+ n-v L+ 2j u dt^, (1) 

and 

A6 = A7rG5a 2 p B , (2) 

where ^f- is the relative temperature variation -with respect to the background 
temperature- along the direction n, and the integral is to be computed from emission 
(e) to observation (o) along the background null geodesies. Initially, this approach 
was designed to study the flat case (scenario I); nevertheless, the potential approxi- 
mation can be also used in the open case (scenario II) provided that we are concerned 
with structures smaller than the horizon scale (Sanz et al. 1996). The first, second, 
and third terms of this equation give the Sachs- Wolfe, the Doppler and the ISW 
anisotropics, respectively. Hereafter, we write A s , A° , or A' to indicate that the 
quantity A has been estimated using only the first, the second or the third term, 
respectively. 

In the linear pressureless approach, the density contrast evolves as follows: 

5(x i ,t) = £^6(x i ,t ), (3) 

where Di(a) describes the evolution of the growing mode of the density contrast 
(Peebles, 1980). The form of the function Di(a) is different in models I and II. 
Hereafter, D lq (a) stands for the functions Di(a) corresponding to our two scenarios. 
The same notation based on the subscript q is also used for other quantities. This 



subscript q is only used along the text to distiguish the model I (q=I) from the model 
II (q=II). Functions Di q (a) and all the quantities having the subscript q are written 
in Appendix A for q=I and q=II. 

The peculiar gravitational potential can be written as follows 

^ x \t) = $(x 1 )^ , (4) 
a 

where function satisfies the equation A$ = B q 5(x l ,to). 

Let us focus our attention on the angular power spectrum of the CMB; namely, 
on the coefficients 

i m=l 

C l=W~ JEW)- (5) 
m=—l 

We begin with the contribution -to these coefficients- of the third term of the right 
hand side of Eq. ([[]), which corresponds to the ISW effect. An appropriate formula 
giving this contribution to the Cgs has been derived. The most useful feature of this 
formula is that, given two k- values and two redshifts, it allows us to obtain a good 
measure of the ISW effect produced by the density perturbation located between the 
chosen redshifts and having scales between the chosen ones. A few words about the 
derivation -similar to the usual derivation of the Ce coefficients of the Sachs- Wolfe 
effect- and characteristics of this formula are worthwhile. 

Since the angular brackets in Eq. (|D stand for a mean performed from many 
realizations of the microwave sky, quantity |a^ m | 2 is first computed for an observer 
having comoving coordinates x p in a reference system attached to the Local Group 
(origin of spatial coordinates) and, then, an average over position x l p is done to get 
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the Ce quantities. 

The equations of the null geodesies passing by the origin of spatial coordinates 



are 



x l = X q (a)e l 



(6) 



Furthermore, in flat cases (as model I), the null geodecis passing by point x l are: 



x l = x l p + \ q {a)e % 



(7) 



while in model II, this last equation is also valid when point P is well inside a sphere 
centred at the Local Group and having the size of the curvature scale. This is because 
-in such a case- the spatial curvature can be neglected. 

Using the third term of the right hand side of Eq. (||), and Eqs. (^), ([|) and (|7|), 
some Fourier expansions lead to the following relation 
AT 



T 



D lg (a) 



da 



where the components of x p are x l , the observer located at P estimates in the 
direction n, and point Q is the intersection between the last scattering surface of P 
and the null geodesic determined by n. 

Finally, from Eq. (|) plus the usual expansions in spherical harmonics and, after 
performing the average in x % , the following angular power spectrum arises: 



a 



^l(k)dk 



(9) 



where 



' je[\ q {a)ka ) — 



D lq {a) 



da 



(10) 



Function P(k) = (\Sk\ 2 ) is the power spectrum of the energy density fluctuations, je 
is the spherical Bessel function of order £, and coefficients T q are given in Appendix 
A for models I and II (all the T coefficients appearing below are also listed in the 
same appendix). 

A similar computation leads to the Cg coefficients corresponding to the first and 
second terms of Eq. (P, which are usually referred to as Sachs- Wolfe and Doppler 
terms. In the Sachs- Wolfe case, these coefficients can be written as follows: 

C l = Kj ^fAK(a)ka }dk (II) 

and, the coefficients of the Doppler term are 

C£ = r; J P(k)j' e 2 [X q (a)ka }dk , (12) 

where j' e (x) = (d/dx)je(x). 

Eqs. (0) - (|T2|) are written in a form which is adequate to perform our nu- 
merical estimates. We define the functions ^\ q {k) = T q k~ 2 P{k)^f q {k) ) ^ q {k) = 
T' 5 q k~ 2 P(k)j 2 [\ q (a)ka ], and ^ q {k) = T q P(k)j' 2 [X q (a)ka }, whose integrals in the 
variable k give , C^ g , and C° , respectively. These definitions will be useful below. 

If the three terms of the right hand side of Eq. (|l|) are simultaneously taken 
into account in order to get (|a^ m | 2 ), the resulting Cg quantities include three crossed 
contributions mixing the ISW, SW, and Doppler effects. We have not found fully 
convincing arguments to neglect these contributions in all the cases and, consequently, 
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they have been systematically estimated using the following formulae: 

2 dk, (13) 
dk , (14) 

2 

dk . (15) 

Since the late ISW effect is produced by inhomogeneities evolving after decoupling, 
quantities C\ q can be estimated using the above pressureless approach for the sources; 
however, the Sachs-Wolfe and Doppler effects are produced by other inhomogeneities, 
which evolve in the recombination-decoupling period and, consequently, a certain 
radiation pressure is acting on the subdominant baryonic component. Taking into 
account that the importance of pressure effects increases as £ does, we only apply our 
pressureles approach to calculate the Sachs- Wolfe, Doppler and crossed coefficients 
in the case £ = 10. This calculation is performed with the essential aim of obtaining 
a rough estimate of the unknown crossed terms and, for this purpose, our approach 
suffices. 

In the £ interval [10,40], we can only expect significant contributions to the Ci 
coefficients coming from: (1) a possible background of primordial gravitational waves 
(this contribution would be almost independent on £ in the interval under considera- 
tion and it is not studied here), (2) each of the three effects considered above and, (3) 
some crossed terms. Other effects as Sunyaev-Zel'dovich, lens anisotropy, nonlinear 
gravitational effects et cettera are not expected to be relevant for these angular scales, 
but for much smaller ones. 

10 



c tq = - 2 / [ve q (k)Ve q (k) 

S I 1 l 

Ve q (k)v eg (k) 



Ctq — 2 



C Tq = -2 J [^q{k)^ q {k) 



The sources of the term C\ have been identified in both scale and position using 
the following definitions: 

-Diq{Z m i n , Z max ) = J — dq{^7 Z m i n , Z max )dk , (16) 



where 



Ciqi^j ^mini ^raax) = a / 

J Zrnin 



d 

dZ 



«o 



dZ . 
(17) 

For Z min = and Z max = Z L , where Z L is the redshift of the last scattering sur- 
face, functions Q q (k, Z min , Z max ) and D\ q {Z min , Z max ) are identical to & q {k) and C' iq , 
respectively. Quantity D' iq (Z m i n , Z max ) can be considered as a measure of the con- 
tribution -to the ISW effect- of the inhomogeneities lying between redshifts Z min 
and Z max ; nevertheless, it is worthwhile to emphasize that the right hand side of 
Eq. fll6|) involves the function (f q (k, Z min , Z max ), which implies that the quanti- 
ties C\ q = D\ q {Q,Z L ) are not the linear superposition of quantities of the form 
D' e (Z min , Z max ), even if these quantities are calculated in disjoint redshift intervals 
covering the total interval (0, Z L ). From Eq. flTBj ) it follows that the contribution 
of each scale to D iq (Z min , Z max ) is measured by the function v\ q {k, Z min , Z max ) = 
T\ q P(k)Q q {k, Z min , Z max )/k 2 . This function measures the contribution of the scale 
k -for the inhomogeneities placed between redshift Z m i n and Z max - to the late ISW 
effect . For Z min = and Z max = Z L , function v iq (k, Z min , Z max ) is identical to func- 
tion /ig q (k) and it weights the contribution of each scale -whatever the inhomogeneity 
location may be- to the ISW angular power spectrum. 
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Our calculations require a value of Z L . Since the Sachs- Wolfe and Doppler effects 
are produced by inhomogeneities located very near the last scattering surface, esti- 
mates of C S iq and C° q based on Eqs. fll^|) and ( |TTD are sensitive to the value of Z L \ 
however, the late ISW effect is produced by inhomogeneities located far from this 
surface (see Section 3) and, consequently, it is almost independent on the assumed 
value of Z L . In order to do the best estimate of the Doppler and SW effects for I = 10 
-allowed by our formalism- we have taken Z L = 1140, which is the redshift corre- 
sponding to ilb = 0.05 and ild = 0.25 according to the formula Z L ~ 1100(r2 rn ,/Sl i) ) 018 
(see Kolb & Turner 1994). Fortunately, we are focusing our attention on the late ISW 
effect, which is almost independent on the choice of Z L . 

3 RESULTS 

Assuming the normalization <t 8 = 1, quantities C\ , C'i 0q , C° , , C^L, C°q , and 
Cio q have been computed in models I (q=I) and II (q=II). For this first normalization, 
Quantity [£(£+ 1)C^ „/27r] 1//2 is shown in the left panel of Fig. 1 (for models I and II). 
The entries 1 and 2 of Table 1 gives [IIOCio^tt] 1 / 2 for all the C\$ quantities. In this 
Table we see that: for I = 10 and model I, the ISW effect is smaller than the Doppler 
and SW ones, while for I = 10 and model II, the ISW and the SW effects are similar. 
A different normalization facilitates some comparisons of the anisotropies appearing in 
models I and II. We have observed that most theoretical predictions based on COBE 
normalization give [£(£ + 1)CV] 1//2 ~ 28 fxK for £ = 10, with a small dispersion around 
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28 \iK. This is true for a wide range of variation of the cosmological parameters: fio? 
fib et cettera. This condition is also compatible with all the observational evidences 
(FIRS, TENERIFE). By these reasons, the ISW effects corresponding to models I 
and II with the normalization [£(£ + l)^] 1 / 2 = 28 fiK for £ = 10 are represented 
in the right panel of Fig. 1, where we see that: (i) the ISW effect corresponding to 
model I (with cosmological constant) is much smaller than that of the model II (very 
open universe), and (ii) the ISW effect of model I is small but it is not negligible. 
Entries 3 and 4 of Table 1 correspond to the second normalization, for which, the bias 
parameter of model I (II) appears to be 0.93 (1.98). This means that the currently 
preferred model (with cosmological constant) leads to a very natural compatibility 
between the CMB observational data and the value a$ ~ 1 extracted from the analysis 
of galaxy surveys. 

In scenario I, the greatest £ = 10 crossed term is the SW-Doppler one (Cf °), 
which is shown in the entries 1 and 3 of Table 1. The remaining crossed terms (SW- 
ISW and Doppler-ISW) are not given because they have appeared to be negligible. In 

SI 

entries 3 and 4 of Table 1, we give the SW-ISW crossed term (C WII ) of the scenario 
II, which is not negligible; however, the terms Doppler-SW and Doppler-ISW can be 
neglected. As it follows from Eqs. (|T3"D - ([15]), any crossed term is proportional to an 
integral (in the variable k), and the function to be integrated can be written as the 
product of two k functions corresponding to the mixed effects. For example, in the 



1/2 



SW-Doppler (SW-ISW) case, we must integrate the product jJ>e„(k) /%,(^) 



1/2 
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( ^i q {k) Heqik) )• In Fig- 2, we display the functions to be multiplied to get 
the SW-Doppler crossed term of model I (left panel) and the SW-ISW term of model 
II (right panel). We have assumed £ = 10 in both models, evidently, these crossed 
terms are not negligible as a result of the existence of a wide enough k interval where 
the positive functions to be multiplied take on large enough values simultaneously. 

Where are located the inhomogeneities producing the ISW effect? This question 
can be answered using Eqs. ([Uj) and flTTD to calculate D iq (0, Z max ) for appropriate 
values of Z max . Results are shown in Fig. 3, where D\ (0, Z max ) is represented as 
a function of Z max in models I (top) and II (bottom). The points of the horizontal 
straight lines of Fig. 3 have the ordinate C\ = -D^(0, 1140) and, consequently, the 
curves D eq (0, Z max ) must tend to the horizontal lines as Z max tends to 1140. In 
the top panel, we see that -0^(0, Z max ) approaches the horizontal lines very quickly. 
Quantity D iT (0,2) is very similar to quantity C\ = -D^(0, 1140), which means that 
the most part of the late ISW is produced by inhomogeneities located at very low 
redshift (Z < 2). This is because the cosmological constant is known to be significant 
only at very low redshifts; before, the universe can be considered as a flat one with 
a negligible cosmological constant, and no ISW effect is expected in this situation. 
In model II (bottom panel), we see that -D^ 7/ (0, Z max ) approaches the corresponding 
horizontal line more slowly than in model I. The most important part of the ISW effect 
is produced by inhomogeneities located at redshift Z < 10, but inhomogeneities at 
Z > 10 also produce a small but appreciable effect. 
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Now, let us look for the spatial scales contributing significantly to the ISW effect 
in the I interval [10,40]. As stated before, the contribution of the scale k to the 
ISW effect -for arbitrary location of the inhomogeneities- is weighted by the function 
ix\ (k). In Fig. 4, this function is represented with solid lines in cases I (top) and II 
(bottom). Left (right) panels correspond to £ — 10 (I — 40). Taking into account 
that, for h = 0.65, the spatial size (diameter in the spherical case) of the structures 
associated to the wavenumber k is \h~ x Mpc, Fig. 4 can be easily interpreted. In 
each panel, the solid lines show a k value for which the function fijJk) reaches a 
maximum. The spatial scale corresponding to this k value is hereafter denoted D*. 
It is the most significant scale for ISW anisotropy production. Solid lines also show 
the existence of a minimum k value where n\ q {k) starts to increase from negligible 
values. The spatial scale associated to the minimum will be denoted D max . Scales 
larger than this maximum one do not contribute to the ISW effect significantly. The 
scales D* and D max corresponding to the four solid lines of Fig. 4 are given in Table 
2. The meaning of these scales is discussed in next section. 

Finally, among all the inhomogeneities located between redshifts Z min = and 
Z ma x, which of them are contributing to the ISW effect? Which are the spatial scales 
of these inhomogeneities? In order to answer this question we have put Z min = 
and various values of Z max into Eq. (|17D ; thus, we have found various functions 
v lq [k, 0, Z max ) of the variable k. Only two of these functions are displayed in each 
panel of Fig. 4. The function corresponding to Z max = 1140 is identical to n\ q {k) 
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and it is drawn with solid lines, while the dotted lines correspond to Z max = 0.5 in 
model I (top panels) and to Z max = 2 in model II (bottom panels). In Fig. 3, we 
see that for these Z max values, a significant part of the ISW effect is produced by 
inhomogeneities located at Z < Z max . The dotted lines of Fig. 4 lead -as the solid 
lines- to new values of D* and D max which are presented in Table 2 and interpreted 
below. These scales correspond to inhomogeneities located at low redshifts smaller 
than the chosen values of Z max 

4 DISCUSSION AND CONCLUSIONS 

In this paper, we have accurately estimated the scales and locations of the inhomo- 
geneities contributing to the late ISW effect. The chosen formalism has facilitated 
our analysis. Now, let us focus our attention on the meaning of the resulting scales 
(see Table 2). They are not the scales of the inhomogeneities (density contrasts) 
producing the effect. According to Eq. ([| ), the ISW effect is produced by the scales 
contributing significantly to the partial time derivative of the peculiar gravitational 
potential and, in the linear regime under consideration, these scales are identical to 
those of the potential itself (see Eq. 

The significant spatial scales of an overdensity and those of its peculiar gravita- 
tional potential are different, this is proved by the relation (fit oc between the 
Fourier transforms of the density contrast 5^ and the peculiar potential (fir. The fac- 
tor 1/k 2 implies that the regions where the potential is significant are more extended 
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than those where the density contrast is not negligible. What is the size of the regions 
where the potential is contributing to the late ISW effect? In order to answer this 



where M(r) is the total mass inside a sphere of radius r. This relation allows us 
to get various important features of the peculiar gravitational potential generated 
by a compensated overdensity. In fact, according to the cosmological principle, any 
overdensity must be compensated at some distance, r c , from its center, at which 
the total mass excess M(r c ) vanishes. This excess also vanishes for r > r c . Then, 
according to Eq. (|18|), the derivative §^r(V c ) vanishes for r > r c and, consequently, the 
potential reaches a minimum constant value, which should be zero to achieve good 
boundary conditions at infinity. We see that the potential of a compensated structure 
tends to zero as r tends to r c ; hence, all the shells forming a certain structure -up to 
compensation radius- would contribute to the ISW effect, although this contribution 
would be small for r values close to r c ; hence, the late ISW effect produced by a 
given structure depends on the way in which it is compensated; namely, it depends 
on the size ~ 2r c of the region where the potential is contributing to the ISW effect. 
This region is hereafter referred as to the pot-region associated to the inhomogeneity. 
Since we are considering linear scales where the peculiar velocities are proportional 
to the gradients of the peculiar gravitational potential, these velocities also vanish for 



question let us consider a spherically symmetric overdensity. In such a case, Eq. (0) 



leads to the relation 
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r > r c and, consequently, the pot-region contributing to the ISW effect is that where 
the peculiar velocities (potential gradients) are significant. A given overdensity would 
contribute to a certain Ce coefficient, if the angular scale subtended by its pot-region 
(not the angular scale subtended by itself) is appropriate. 

The compensation of cosmological objects is a statistical phenomenon and, con- 
sequently, structures of the same type (for example various GAL structures) could 
be compensated at different distances from their cores; therefore, although the com- 
pensation radius of the Great Attractor has been estimated to be r° c A ~ lOO/i -1 Mpc 
(study of the velocity field around the GA), other GAL structures could compensate 
at other distances, perhaps at distances of a few times r^ A . Voids and Abell clusters 
would compensate at distances of a few tens of Mpc from the central region. Taking 
into account these considerations we are going to interpret the results summarized in 
Table 2. 

For model I and £ = 40 (£ = 10), pot-regions with radius larger than 250k' 1 Mpc 
(1000/i -1 Mpc) are not contributing to the late ISW effect. Those having radius 
larger than 50/i _1 Mpc (200/2T 1 Mpc) and located at Z < 0.5 do not contribute ei- 
ther. The maximum effect is produced by pot-regions with radius close to 100/iT 1 Mpc 
(300/ir 1 Mpc) and located between redshifts 0.5 and 2. and, finally, the maximum ef- 
fect produced at Z < 0.5 is due to pot-regions having about 40/i _1 Mpc (140/i _1 Mpc) 
radius. This means that, in model I, GAL structures with r c ~ 100/i" 1 Mpc produce 
the maximum contribution to C40 (the smallest of the C' n coefficients, see the top 
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panel of Fig. 1). GAL objects with sizes r c ~ 140/i _1 would contribute to Cio when 
located at Z < 0.5 and GAL objects with pot-regions of a few times lOO/i -1 Mpc 
would contribute to all the C' n coefficients from I = 10 to I = 40. Pot-regions with 
radius of ~ 40/i _1 Mpc and located at Z < 0.5 contribute to C± . 

In model II, a similar study has been developed. For £ = 40 (£ = 10), the following 
conclusions can be obtained: (i) for arbitrary locations, the most large compensation 
radius contributing to the late ISW is 650/i _1 Mpc (2500/i -1 Mpc), (ii) for structures 
located at Z < 2, the most large radius is 125/i _1 Mpc (500/iT 1 Mpc), (iii) for arbitrary 
locations, the most large contributions to the late ISW come from compensation 
radius of 300/i _1 Mpc (1000/i" 1 Mpc), and (iv) for structures located at Z < 2, the 
most large contributions correspond to radius of ~ llO/i -1 Mpc (~ 350/1^ Mpc). In 
model II, the scales are larger than those of model I. Scales of a few times ~ 10/z 1 Mpc 
do not play any role. GAL structures with r c ~ lOO/i^ 1 Mpc only contribute to C 40 
if located at Z < 2. GAL objects compensated at radius around SOOh' 1 Mpc would 
play an important role in generating the C' n coefficients, in particular, for i = 40. 
This conclusion is in agreement with previous claims about the possible relevance of 
GAL structures in generating the late ISW effect in open universes (Arnau, Fullana 
& Saez, 1994; Saez, Arnau & Fullana, 1995). The GAL objects simulated in those 
studies (based on TB) undergo effective compensations at distances of a few hundred 
of Megaparsec from their cores and, in agreement with the results of this paper (Table 
2), this type of structures would be contributing significantly to the late ISW effect. 
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It is worthwhile to emphasize that we have discussed the contribution to the late 
ISW effect of great cosmological structures (which produce peculiar velocities up to 
distances of tens or hundreds of Mpc). A linear approach suffices to estimate the 
potential (also the peculiar velocities) produced by these structures (GAL objects, 
voids et cettera). We have never considered the Rees-Sciama effect produced by 
strongly nonlinear substructures lying inside the Great Attractor and other extended 
inhomogeneities. Such an effect would produce CMB anisotropy on smaller angular 
scales and its estimate would require other nonlinear approaches. 



APPENDIX A 



Some quantities used in this paper have the subscript "q" . Here, the explicit form of 
these quantites is given for q = I (model I of Section 1) and q = II (model II). We 
summarize the information as follows: 
MODEL I (q = I) 

The growing mode of the scalar energy density fluctuations is 



D u (a) 



x 



I 1/2 



+ x z 



X 



f 

Jo 



-3/2 



dy , 
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where 



x 



2Q> 



1/3 



[i + z]- 1 



(20) 



The constant B q is 



(21) 
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The function X q (a) can be written as follows: 



A/ (a) = K(a)H 1 , 



(22) 



where 



K(a) 



db 



(23) 



L (tt m0 b + n x b*y/* ■ 

I S D SD SI DI 

Now, we give the coefficients Tj, Fj, Fj , Fj , T 7 , and Tj defined in Section 2. 



n 



mO 



Di(0) 



(24) 



2ir 



D 1 (L)n m0 (l + Z L ) 
D 1 (0) 



(25) 



D _ 2H 

1 ~ ~v~ 



Di(L) 
DM 



n m0 (i + z L ) . 



(26) 



In order to derive these formulae, the cosmological constant has been assumed to 
be negligible at Z L (see comments of Section 3) and, consequently, as Q m o tends to 

S D 

unity, quantities V j and T 7 tend to the right values corresponding to a flat universe 
without cosmological constant (Hq/2h and 2H^/n(l + Z L )). 
MODEL 11 (q = II) 

The same quantities as in model I are now listed: 



where 



DuM = i + - ( + 3(1 c + 3/ f /2 /n[(i + C ) 1/2 - C 1/2 ] , 
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(27) 



(28) 



B 



11 



— O 
2 



# £iK)(l-fio) 3/2 



-i 



(29) 



X n (a) = 2tanh(y/2) 



(30) 



where 



y = cosh 1 



2-n c 
n 



cosh 1 



2(1 - n f/ 2 H 

fin 



a + 1) 



(31) 



1 // 



tt(1 - fi )^f (ao) 



(32) 



n 



2tt 



°nl(i + z L f 



DIM 
D?(oo) 



(33) 



2fi 



'd£>i' 



(34) 



" 7T(l-fio)A 2 (0)(l + ^) V rf « / L ' 

As fio tends to unity, function Di(a) tends to the growing mode of a flat back- 
ground, which is proportional to a. Taking into account this fact and Eqs. (|33"|) and 
(p4j) , one easily concludes that -as the universe approaches a flat one quantities T^-j 
and T^j tend to Hq/2tt and 2Hq/tt(1 + Z L ), respectively. These limit values coin- 

S D 

cide with the well known values of r 7J and T /7 corresponding to the flat background 
without cosmological constant. 
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FIGURE CAPTIONS 



FIG. 1.— Each panel shows the quantity [£(£+ l)Cj, q /2ir} 1 / 2 (in /iK) as a function of 
I. Left (right) panel corresponds to the normalization a% — 1 ([£(£ + l)Cy 1//2 = 28 /z-fT 
for £ = 10). 

FIG. 2.— Left panel shows the quantity /4oj x 10 5 defined in the text as a function of A; 
for the SW (solid) and Doppler (dotted) effects. In the right panel, quantity /ijg/ x 10 4 
is plotted vs. k for the SW (solid) and ISW (dotted) effects. These functions must be 
multiplied to get the corresponding crossed contributions to the CMB angular power 
spectrum. 

FIG. 3.— Top panel shows the quantity D' £I (0, Z max ) defined in the text as a function 
of Z max for I = 10 (solid, not horizontal) and £ = 40 (dotted, not horizontal). These 
curves approach the horizontal lines of the same type as Z max increases, crossing them 
at Z max = 1140. Bottom panel has the same structure, but it shows the quantity 
D UI (0, Z max ). 

FIG. 4.- Top left: plot of v WI {k, 0, Z max ) (see text) vs. k for Z max = 1140 (solid) 
and Z max = 0.5 (dotted); top right: the same as in top left panel for the quantity 
v' 40I (k, 0, Z max ); bottom left: plot of v 10II (k, 0, Z rnax ) vs. k for Z max = 1140 (solid) and 
Z min = 2 (dotted); and bottom right: the same as in bottom left for v mu {k, 0, Z max ). 
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TABLE 1 

[1(1 + 1)Q/2tt] 1 /2 (in fiK) FOR £ = 10 



MODEL a 8 ISW SW Doppler SW-Doppler SW-ISW Total 

I 1.00 6.75 21.85 12.63 1.25 26.15 

II 1.00 38.18 38.03 12.63 12.30 55.34 

I 1.07 7.23 23.39 13.52 1.34 28.00 

II 0.51 19.32 19.24 6.39 6.22 28.00 
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TABLE 2 

PRESENT SIZES OF THE INHOMOGENEITIES 
PRODUCING THE LATE ISW EFFECT 



MODEL 


i 


'Z'max 


D m ax 

h- 1 Mpc 


D* 

h,- 1 Mpc 


I 


10 


1140 


2000 


600 


I 


10 


0.5 


400 


280 


I 


40 


1140 


500 


200 


I 


40 


0.5 


100 


85 


II 


10 


1140 


5000 


2000 


II 


10 


2. 


1000 


700 


II 


40 


1140 


1300 


600 


II 


40 


2. 


250 


220 
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